A Classification Scheme for Phenomenological Universalities in Growth Problems 
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A classification in universality classes of broad categories of phenomenologies, belonging to dif- 
ferent disciplines, may be very useful for a crossfertilization among them and for the purpose of 
pattern recognition. We present here a simple scheme for the classification of nonlinear growth 
problems. The success of the scheme in predicting and characterizing the well known Gompertz, 
West and logistic models, suggests to us the study of a hitherto unexplored class of nonlinear growth 
problems. 



Or 
I ' 

o ■ 



C/3 

o 



> 

in 
o 

o 

O 

o 

^3 : 

Or 



X 



Present efforts towards the understanding of complex 
systems in physics, biology, economics and social science 
require complementary microscopic and macroscopic de- 
scriptions. In fact, due to the complexity of the underly- 
ing dynamics and the unbounded variety of external con- 
ditions, a fundamental approach is missing. Microscopic 
models depend on such a large number of parameters 
that they often lose almost any predictive power, even 
when the calculations do not become forbiddenly diffi- 
cult or time consuming. On the other hand, macroscopic 
descriptions are often inadequate and do not take advan- 
tage of the enormous progress that has been achieved at 
the microscopic level in recent years. An intermediate 
(mesoscopic) approach 0, 0, ma Y be very fruitful, but 
a bridging among the various levels Q is not always easy 
to accomplish. 

A different approach has consequently emerged for the 
treatment of problems , which do not directly require 
a detailed description of the system to be investigated. 
The idea is to exploit the spectacular advancement of 
interdisciplinary research, which has taken place in the 
last two decades or so, involving e.g. the relevance of 
scale laws, complexity and nonlinearity in virtually all 
disciplines. 

In this context many patterns have been discovered, 
which are remarkably similar, although they concern 
completely different phenomenologies. This is hardly sur- 
prising, since often the "background" mathematics is the 
same. We shall call them "phenomenological universal- 
ities" 0, in the sense that they refer to a "transversal" 
generality (not to a uniformly general behaviour within 
a given class of phenomena) . 

As examples of universality we can quote the "life's 
universal scaling laws" 0, which will be discussed later, 
and the "universality of nonclassical nonlinearity" Q. 
The latter suggests that unexpected effects, such as those 
recently discovered by P. Johnson and collaborators || 
and called by them "Fast Dynamics", may be found as 
well, although possibly with quite different manifesta- 
tions in other fields of research. 



A reliable macroscopic analysis of a complex system 
requires two fundamental ingredients: non linearity and 
stochasticity. Non linearity is more fundamental because 
the stochastic behaviour requires a non linear dynamics. 
Therefore non linearity must be considered as the fun- 
damental feature of these systems and in this letter we 
consider general growth problems based on this crucial 
aspect. We shall show that different "degrees of nonlin- 
earity" (as specified below) correspond to various growth 
patterns, which can be systematically classified. 

For this purpose, let us consider the very broad class 
of growth phenomena, which may be described by the 
simple law: 



dY(t) 
dt 



a(t)Y(t) 



(1) 



where a(t) represents the specific growth rate, which 
may vary with time, of a given variable Y(t). By intro- 
ducing the nondimensional variables t = a(0)t, y(t) = 
Y(t)/Y(0) and a(r) = a(i)/a(0), Eq.(l) becomes: 



dr 



a(r)j/(r) 



(2) 



with y(0) = a(0) = 1. By defining the time variation of 
a(r) through a function $(a): 



*(a) = - 



da(r) 
dr 



(3) 



we obtain a system of two differential equations, which 
may generate a variety of growth patterns, according to 
the explicit form of ^(a) , and is usually analyzed by the 
standard fixed points and characteristic curves methods 

In this contribution we are not directly interested in 
this aspect, but we wish to show, instead, how the non- 
linear terms in $(a) affect the growth dynamics process. 

We assume that |a(r)| < 1 and expand ^(a) in power 
series 



*(a) = £~ =0 &„a n 



(4) 
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in which we retain only a limited number of TV + 1 
terms. Borrowing from the language of phase transitions 
[Io| , we define, as belonging to the phenomenological 
universality class of order N (which we shall call UN, 
N=l,2,...), the ensemble of all the phenomenology, which 
may be suitably described by truncating the series at the 
power n=N. In the following we shall analyse in detail 
the classes Ul, U2 and U3 and provide a description of 
their nonlinear properties. 

The "linear" behaviour of the system corresponds to a 
constant specific growth rate, i.e. $(a) = ( or b n = for 
any n). Then y(r) follows a purely exponential law. Also 
the case bo ^ with all b n = for n > 1, can be easily 
shown to lead to an exponential growth. Since we are 
interested only in the nonlinear effects, we shall assume 
bo = 0. This does not cause any loss of generality, since 
one can always expand $ in the variable (3 = a — c, where 
c is a solution of H%L b n c n = 0. In the (3 expansion the 
coefficient of f3° vanishes. Likewise, again without any 
loss of generality , we can set b\ = 1, as one would have 
from an expansion in the variable 7 = a/61. 

In order to study the various classes of universality 
and obtain the corresponding differential equations and 
solutions, we write from Eqs. (2) and (3): 



$(a) 



da 



ay 



from which it follows: 



Iny 



ada 



const 



(5) 



(6) 



By solving the previous equation with respect to the vari- 
able a(r) and then substituting into Eq. (2), one obtain 
the differential equation characterizing the class. The in- 
tegration constant can be easily obtained from the initial 
conditions. 

Let us then start by considering the class Ul, i.e. with 
N=l. From Eq. (6) and ^(a) = a , it immediately 
follows: 



dy 
dr 



= y- ylny 



with the solution 



y — exp[l — exp (— r) 



(7) 



(8) 



Eq.(7) represents the "canonical" form of Ul differential 
equations and cor resp onds to the Gompertz law, orig- 
inally introduced 1 1 11] in actuarial mathematics to eva- 
lute the mortality tables and, nowdays, largely applied 
to describe economical and biological growth phenomena. 
For example, the Gompertz law gives a very good phe- 
nomenological description of the tumor growth pattern 
E3> [l| an d it can be related to the energetic cellular 
balance |l4j . It is remarkable that it does not contain any 
free parameter ( except for the scale and linear parame- 
ters which have not been included, as discussed before), 



i.e. all Gompertz curves are (under the mentioned pro- 
viso) identical. 

Let us now turn our attention to the class U2, i.e. N=2. 
From Eqs. (6) and (3) and ^(a) = a + ba 2 , where b = 62 
, it follows 



dr 



a 2 y p - fay 



(9) 



where a 2 = (1 + b)/b, p = 1 — b and = 1/6 with the 
solution 



y = [1 + 6 — 6 exp (— r 



,1/6 



(10) 



By identifying y with the mass of a biological system, y = 
to, and defining the asymptotic mass (too = yo = 1) 



M = lirrir^^m^T) = (1 + 6) 



(11) 



it is easy to show that Eqs. (9) and (10) correspond 
to the well known allometric West equation for the case 
p = 3/4 [l5j |. In their ontogenetic growth model, m rep- 
resents the mass of any living organism, from protozoa 
to mammalians (including plants as well). By redefin- 
ing their mass and time variables z = 1 — (y/m) b and 
8 = —t + Inb — blnM they obtain the very elegant pa- 
ramcterless universal law 



z = exp{—t 



(12) 



which fits well the data for a variety of different species, 
ranging from shrimps to hens to cows. It is interesting to 
note that, in a subsequent work 0], West and collabo- 
rators give an interpretation of 9 as the "biological time" 
, based on the organism's internal temperature. 

An extension of West's law to neoplastic growths has 
been recently suggested by C. Guiot, P.P. Delsanto, T.S. 
Deisboeck and collaborators 0, • Although an un- 
ambigous fitting of experimental data is much harder in 
tumors (except for cultures in vitro of multicellular tumor 
spheroids), the extension seems to work well . Of course, 
particularly in vivo, other mechanisms must be taken into 
account, such as the pressure from the surrounding tissue 
[Tij. Another important issue is the actual value of the 
exponent p, which has been the object of a strong debate 
[lo| • Recently C. Guiot et al. have proposed that p 
may vary dinamically with the fractal nature of the input 
channels (e.g. at the onset of angiogenesis). 

Although it is not obvious from a comparison between 
Eq. (7) and Eq. (9), Ul represents a special case (6 = 0) 
of U2, as it obviously follow from the power expansion of 
<!> (which has 6 = in Ul). This can be verified directly 
by carefully performing the limit 6 — * in Eq. (10) . In 
fact it is interesting to plot y vs. r in a sort of phase 
diagram ( see Fig. 1) . 

This leads to a very suggestive interpretation of Eq. 
(9). Having added a term to the $(<z) expansion, we 
gain, in U2, the possibility of adding a "new" ingredi- 
ent, which turns out to be a different dimensionality of 
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FIG. 1: - Growth curves belonging to the class U2. From 
the top to the bottom the values of the parameter b are 
-0.25,-0.1,0.1,0.25,0.5 respectively. The solid curve (b = 
0,p = 1) corresponds to the Gompertzian (Ul), while the 
dashed one refers to the value proposed in Q p — 3/4 
(6=1/4). 



the "energy flux" i.e. input, output and consumption 
(metabolism). E.g. the first term on the RHS of Eq. 
(9) may be related [12] to the premise that the tendency 
of natural selection to optimise energy transport has led 
to the evolution of fractal-like distribution networks with 
an exponent p for their terminal units vs. an exponent 1 
for flux mechanisms related to the total number of cells. 
When b=0, p=l and we lose the new ingredient, thus 
falling back into Ul. 

This is confirmed also by considering the logistic equa- 
tions, corresponding to eq.(9) with negative b. The usual 
logistic equation is obtained for p = 2. As well known in 
population dynamics |2.'"'J| , in this case the new ingredient 
is the competition for resources. 

Finally we consider the class U3. Writing 



$(a) = a(l + ba + ca 2 ) 



from Eq. (6) it follows 



/ 



da 



1 + ba 



K — Iny 



(13) 



(14) 



In this case there are three subclasses, U31, U32 and U33, 
corresponding to A = 4c — b 2 > < . For brevity we limit 
ourselves to report here the canonical equation for U31, 
i.e. when A < 0: 



(It 



a-zy - (i-iV 9 + 73 



dr 



(15) 



where d = y/^E ,p = 1 - d,K = (d - 3c)/ (d + 3c), 
a 3 = (d-c)/2c and /3 3 = K(d+c)/2c and 73 = K/(l-d). 
It is interesting to observe that, in the same way that 
U2 adds ( with respect to Ul) a term with a different 
dimensionality to the energy flux contribution, U3 adds 
such a term (the last one in Eq. 15) to the growth part. 



To conclude, we have developed a simple scheme, 
which allows the classification in nonlinear phenomeno- 
logical universality classes of all the growth problems, 
which can be described by Eqs. (2) and (3). We have 
found that the first class Ul corresponds to the Gompertz 
curve, which has no free parameters (apart from scale 
and linear ones). The second class U2 includes all the 
Westlike and logistic curves and has a free parameter b: 
when b=0 we fall back into Ul (Gompertz). The success 
of the scheme in obtaining the classes Ul and U2 when 
one or two terms are retained in the expansion of &(a) 
has suggested to us to investigate the class U3, which is 
generated by simply adding one more term (see Eq. 13). 
To our knowledge, this class has never been investigated 
before. A remarkable result is that each new class adds a 
new "ingredient" (or growth mechanism). E.g. U2 allows 
for the possible presence of two dimensionalities in the 
energy flux. U3 extends such a possibility to the growth 
term ( the time derivative). In addition to its intrinsic 
elegance [24| the concept of universality classes may be 
useful for several reasons of applicative relevance. In fact 
it greatly facilitates the crossfertilization among different 
fields of research by implicitly suggesting that a method 
of analysis, which is proven advantageous in one study, be 
tried and eventually adopted in others. Also, if an unex- 
pected effect is found experimentally in a field, similar ef- 
fects "mutatis mutandis" should also be sought in similar, 
although unrelated, experiments in other fields. Finally, 
if a detailed study is performed to recognize the patterns 
that are characteristics of the most relevant classes (and 
subclasses) , this could greatly help in classifying and fit- 
ting new sets of experimental data independently of the 
field of application. 
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